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Abstract 

The paper discusses the similarities and the differences in the mathematical 
theories of the steady Boltzmann and incompressible Navier-Stokes equations posed 
in a bounded domain. First we discuss two different scaling limits in which solutions 
of the steady Boltzmann equation have an asymptotic behavior described by the 
steady Navier-Stokes Fourier system. Whether this system includes the viscous 
heating term depends on the ratio of the Froude number to the Mach number of 
the gas flow. While the steady Navier-Stokes equations with smooth divergence-free 
external force always have at least one smooth solutions, the Boltzmann equation 
with the same external force set in the torus, or in a bounded domain with specular 
reflection of gas molecules at the boundary may fail to have any solution, unless the 
force field is identically zero. Viscous heating seems to be of key importance in this 
situation. The nonexistence of any steady solution of the Boltzmann equation in 
this context seems related to the increase of temperature for the evolution problem, 
a phenomenon that we have established with the help of numerical simulations on 
the Boltzmann equation and the BGK model. 


Introduction 

The Boltzmann equation and its hydrodynamic limits have been widely 
studied in the time-dependent regime. The Cauchy problem for the Boltz- 
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mann equation is discussed in [Ml ED 133 US]. Hydrodynamic limits of the 
Boltzmann equation are analyzed by various methods in BD ns El ESI IS 0 
[9l0|35lE6ll23lll2l|2ll[25l|33] — see also |53j for a nice introduction to the 
mathematical analysis of hydrodynamic limits of the Boltzmann equation. 

Some of the results known in the case of the Cauchy problem set on a 
spatial domain which either the Euclidean space or the periodic box 
have been extended to the case of the initial-boundary value problem in a 
domain Q of R^: see jlOl |39l 0 ED — see also |33j for a more synthetic 
presentation of this material as well as further results. 

By comparison, the mathematical literature on the analogous problems 
in the case of steady solutions is much more scarce. A great collection of 
asymptotic and numerical results on the Boltzmann equation in the steady 
regime can be found in the books |45ll46j . As for the mathematical analysis of 
the boundary value problem for the Boltzmann equation, the main references 
are [271128] (see also [26] ) , together with the more modern references [D HD- 
Steady solutions of the Boltzmann equation for a gas flow past an obstacle 
have been investigated in detail in [Hp 152]. 

There are striking analogies between the Boltzmann and the incompress¬ 
ible Navier-Stokes equations in three space dimensions — in the words of 
P.-L. Lions [33] “[...] the global existence result of [renormalized] solutions 
[...] can be seen as the analogue for Boltzmann’s equation to the pioneering 
work on the Navier-Stokes equations by J. Leray”. This analogy is at the 
origin of the program outlined in and carried out in [231 [25] — see 

also |33j for the extension to a more general class of collision kernels. 

It has been known for a long time that the regularity theory of solu¬ 
tions of the incompressible Navier-Stokes equations is much easier in the 
steady than in the time-dependent regime. For instance, in space dimension 
3, steady solutions of the Dirichlet problem for the incompressible Navier- 
Stokes equations have the same regularity as their boundary data and the 
external force field driving them: see Proposition 1.1 and Remark 1.6 in 
chapter II, §1 of [39|- At variance, it is still unknown at the time of this 
writing whether Leray solutions of the Cauchy problem for the Navier-Stokes 
equations in space dimension 3 propagate the regularity of their initial data 
— see Problems A-B in m- 
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One striking difference between the steady and the time-dependent prob¬ 
lems for the incompressible Navier-Stokes equations is the uniqueness theory, 
and its relation to the regularity of solutions. Smooth solutions of the time- 
dependent incompressible Navier-Stokes equations in space dimension 3 are 
known to be uniquely determined by their initial data (and driving force 
field) within the class of Leray weak solutions, a remarkable result proved 
by Leray himself (see §32 in [32]). On the contrary, it can be proved that 
bifurcations do occur on the steady problem for the incompressible Navier- 
Stokes equations, leading to nonuniqueness results. Such a nonuniqueness 
result for the Taylor-Couette problem has been proved in (5^ — see also 
Chapter II, §4 in |49j . Analogous bifurcations for the Boltzmann equation 
have been observed numerically in j48| . and mathematically in [2]. 

These considerations suggest studying the mathematical theory of exis¬ 
tence and regularity for steady solutions of the Boltzmann equation driven 
by an external force field. In particular, does the assumption of a steady 
regime simplify regularity issues, as in the case of the Navier-Stokes equa¬ 
tions? 

As we shall see, there are strong similarities between the Boltzmann and 
the Navier-Stokes equations in the steady regime. After reviewing the basic 
structure of the Boltzmann equation in section 1, we propose in section 2 a 
formal derivation of two variants of the incompressible Navier-Stokes-Fourier 
system from the Boltzmann equation under appropriate scaling assumptions. 
Section 3 discusses the steady Navier-Stokes and Boltzmann equation in the 
periodic setting, which leads to a striking difference between both models. 
Section 4 provides a physical explanation for this difference, based on numer- 
icla simulations of the evolution problem. After a brief section 5 summarizing 
our conclusions, some computations involving Gaussian averages of certain 
vector and tensor fields have been put together in an appendix. 

Professor Tai-Ping Liu is at the origin of some of the most striking 
results on the mathematical analysis of the equations fluid dynamics — his 
work |37| on the compressible Euler system, and his analysis of the stability 
of the Boltzmann shock profile in collaboration with S.-H. Yu [38| have had 
a lasting impact on our field. We are pleased to offer him this modest 
contribution on the occasion of his 70th birthday. 
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1. The Steady Boltzmann Equation with External Force Field 


The Boltzmann equation with external force field / = f{x) G is 
posed in a bounded, convex spatial domain O C with smooth boundary 
dQ. The outward unit normal vector field on dO, is denoted by n^- According 
to §1.9 in |46] . its dimensionless form is 

Ma^ 1 q 

V ■ VxF + -z^f ■ VyF = :^C{F) , X G n , u G R^ , (1.1) 

rv 


where 


CiF) := 


R3xS2 


(f'f: 


FF^:)\{v — u*) • uj\dv^duj , 


with the usual notation 


F = F{x,v), F^ = F{x,v^), F' = F{x,v'), 


and F^ = F{x,vJ , 


assuming that 

v' = V — (v — v^) ■ ojoj , and = u* + (u — u*) • uu . 


The dimensionless numbers Ma, Fr and Ki are respectively the Mach, Froude 
and Knudsen numbers. We recall the definitions of these dimensionless num¬ 
bers: 


Ma = 


Up 


Fr = 


Up 

VFpU' 



where R is the specific gas constant, henceforth set to one for simplicity. In 
these formulas, Uq, Tq, Fq, Lq are respectively the reference speed, tempera¬ 
ture and external force in the gas, while Lq is the reference length scale and 
£q is the mean free path of the gas molecules at the reference state. 


We recall that the local conservation laws of mass, momentum and en¬ 
ergy for the collision integral are 



C{F)dv 



ViC{F)dv 



v\^C{F)dv = 0 


a.e. on 


for i = 1,2,3, provided that E > 0 is measurable and decays rapidly enough 
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as |u| —)• cx) — so that, for instance 



(1 + |n|^)(l + \ f\)Fdvdx < oo . 


(See for instance chapter I, §4, Corollary 1 in [10].) This implies the local 
mass, momentum and energy balance identities, which are the differential 
identities 


diVa; / vFdv = 0 , 
4r3 


div^, 

diVa; 


f ^ , Ma 
/ n (8) vFdv = —^ 
/r 3 Fr2 


/R3 


v^\v\‘^Fdv = 


Ma^ 


/ / Fdv, 

Jr3 

f ■ [ vFdv . 
J-R3 


mass 


(momentum) 

(energy) 


are satisfied by all the solutions of the Boltzmann equation having the decay 
property mentioned above. Notice that this property implies that 


[ [ F\f\ds{v)dx^0 
JfiJ |ti|=_R„ 

for some sequence Rn —>• oo. (The notation ds designates the surface element 
on the sphere of radius Rn centered at the origin.) 

Integrating fnrther in x and applying Green’s formula leads to the iden¬ 
tities: 





tan Jr? 



Fv ■ nxdvdS{x) = 0 , 


R3 

vF{v ■ nx)dvdS{x) 
pF(u • nx)dvdS{x) 


Ma^ 

Ma^ 

FV^ 



fFdxdv , 

V ■ fFdxdv , 


(mass) 

(momentum) 

(energy) 


where dS is the surface element on did, which are the global balance laws of 
mass, momentum and energy. 

Multiplying both sides of the Boltzmann equation by InT + 1 leads to 
the identity 


V ■ Vx{FlnF) + 


Ma^ 

FY^ 


/•V^(FlnF) 


1 


C{F){lnF + l). 
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Integrating in v leads to the local form of Boltzmann’s H theorem (see for 
instance chapter I, §4, Corollary 1 in |20j l 

diva; [ vF\nFdv = — [ C{F) In Fdv < 0 , 

Jr3 

assuming again that F decays rapidly enough as |u| —)• oo — for instance 


OxRS 


(1 + |u|)(l + \ f\)FlnFdvdx < oo 


so that 

Rn [ [ F\\nF\\f\ds{v)dx^O 

Jo. J |u|=_R„ 

for some sequence Rn oo. 

Integrating further in x, one obtains the global form of Boltzmann’s H 
theorem 


[ [ F In Fv ■ nxdvdS{x) = 

Idfi Jr^ Kr 


OxR3 


C{F) In Fdvdx < 0 , 


Boltzmann’s H theorem asserts that the inequality above is an equality if 
and only if F is a local Maxwellian, i.e. is of the form 


F(x, u) 1 


with the notation 

P \v-u\^ 

= (27r0)3/2^ • 

(See for instance chapter I, §§5 and 7, Corollary 2 in [lUj-) ■ 


( 1 . 2 ) 


2. The Navier-Stokes Limit for the Boltzmann Equation 

In this section, we explain how two variants of the Navier-Stokes-Fourier 
system can be derived from the Boltzmann equation. The exposition is for¬ 
mal and follows the style adopted in m s]- The difference between the 
limiting systems comes from the different scaling assumptions on the Froude 
number. Specifically, we are concerned with those variants of the incom¬ 
pressible Navier-Stokes-Fourier system which do, or do not include the vis- 
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cous heating term. This particular feature of the fluid dynamic limit of the 
Boltzmann equation is of key importance for the discussion in the present 
work. 


2.1. Prom Boltzmann to Navier-Stokes-Fourier 


In this section, we assume the following scaling, where e > 0 is a small 
parameter: 

Ma = Ki = e , andFr = \/e. (2-1) 

Hence, the scaled Boltzmann equation takes the form 

V ■ VxFe + ef ■ VvFe = -C{F^) , X G H , u G . 

e 

In addition, write the Helmholtz decomposition of the external force field as 

/(x) = -V$(x) + efs{x ), with div fs = 0. 

(In other words, we assume that the divergence free component of the exter¬ 
nal force is small compared to its curl free component.) With this additional 
assumption, the scaled Boltzmann equation becomes 

V ■ VxFe — eV<I>(x) • + e^fs{x) ■ = -C{F ^), x G H , x G . 

(2.2) 

Seek Fe in the form 


Fe{x,v) = M{v) + egeix,v)^ , 


assuming that 


with the notation 



QeMdxdv = 0 , 


M AI(i^o,i) ’ 


and 


1 


in 


-x. 


(2.3) 


(2.4) 


(2.5) 
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In other words, the distribution function is sought in the form of a perturba¬ 
tion of the order of the Mach number Ma = e about the uniform Maxwellian 
M. Except for the scaling assumption on the external force, this is exactly 
the same scaling assumption as in las!. 


In terms of g^, the scaled Boltzmann equation (2.2) takes the form 


ev • • V,{Mg,) + e^M-^f,{x) • + e^,)) 

( 2 . 6 ) 

where 


Cg := -M-^DC{M) ■ {Mg ), Q{g, g) := M-^C{Mg ), (2.7) 


are respectively the linearized collision integral at M and the collision inte¬ 
gral intertwined with the multiplication by M. (The notation DC{M) ■ {Mg) 
designates the differential of the collision integral C evaluated at M, and ap¬ 
plied to the variation Mg of distribution function.) The quadratic operator 
Q defines a unique bilinear symmetric operator, also denoted Q, by the 
polarization formula 

Qif,g) ■=\{Q{f + 9,f + 9)- Q{f, f) - Qig,g)) ■ 


In other words, 


^9= {g + g^ - g'- gi)\{v - v^) ■ uj\M^dv^duj, 


R3xS2 


while 


Q{f,9) = h 


2 I I ifg'* + f'*9 - fg* - f*9)\{v - v^) ■ uj\M^dv^duj. 
R3xS2 


Henceforth, the integration with respect to the Gaussian weight M is 


denoted as follows: 


((/>):= f c{){v)M{v)dv. 
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Theorem 1. Let be a family of solutions of the scaled Boltzmann equation 


(2.2), whose relative fluctuation g,. defined in (2.3) satisfies 


g 

weakly in Lf{Ll x R^; (1 + \v\‘^)Mdvdx), and 

Q{9e,9e) Qig,9) 

weakly in L^(R^; (1 + \v\'^)Mdv)), while 

eigeVfiv)) 0 

weakly in L‘^{VL) for each (f G L^(R^; (1 + \v\‘^)Mdv). 
Then 

1 n „,|2 


g{x, v) ■.= 6 + u{x) ■ V + 0{x)^{\vf — 5) 


where 


and {u, 9) is a solution of the incompressible Navier-Stokes-Fourier system 
div w = 0 , 


div(n®^) + Vp = izAu + + fs , 


_ I div(u0) = kA6 — tt • Vd>. 


( 2 . 8 ) 


The values of the viscosity v and heat diffusivity k are determined implicitly 


in terms of the collision integral, by formulas (2.15) and (2.16) 


We recall the following fundamental result. 

Lemma 1. The operator L is an unbounded self-adjoint operator on the 
Hilbert space L‘^(R?-,Mdv) with domain L^(R^; (1 + \v(^)Mdv). Moreover 

£ > 0 and Ker£ = span{l, vi, V 2 , V 3 , |u|^} . 
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Finally, C satisfies the Fredholm alternative: one has 

Ran£ = (Ker£)-*- . 

This has been proved by Hilbert in 1912 (see for instance [20], chapter 
III, §§4-5 and [TT] . chapter IV, §6). 


Proof. The proof of Theorem is rather involved; we follow the discussion 

in [Sill]. 

Step 1: asymptotie form of the number density fluctuations 

Assuming that g,. ^ g weakly in x Mdxdv), we pass to the 


limit in the sense of distributions in both sides of the equality (2.6), and 
obtain 

Cg = 0. 

Thus g is of the form 

g{x, v) = p{x) -b u{x) ■ V + 6{x)^{\v\'^ — 3). (2-9) 


Step 2: divergence-free and hydrostatic conditions 

Multiplying both sides of the scaled Boltzmann equation by and 
integrating in v leads to 

dwa:{vge) = e [ divy{M{g^'^^{x) - ege)fs{x)))dv = 0. 

JR3 

Passing to the limit in both sides of this equality in the sense of distributions 
as e —>■ 0, we get 

diVa; u = divxivg) = 0 . (2-10) 

Mutiplying both sides of the scaled Boltzmann equation by ^Mv and 
integrating in v leads to 

divx{v®‘^ge) = -eV^{x){gfi + efs{x){{Z,,e^'^‘'^'^ + egf)) • 
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Passing to the limit 


Passing to the limit 


as e —>• 0 shows that 


Vx(/0 + 0) = divxiv^'^g) = 0 . 


in (2.4) shows that 


p{x)dx = / {g){x)dx = 0 , 

1 Jq. 


so that 


p{x) + 9{x) = 9 , 


with 9 := 


1 

W\ 



9{x)dx . 


( 2 . 11 ) 


( 2 . 12 ) 


Step 3: motion equation 

Mutiplying both sides of the scaled Boltzmann equation hy -^Mv and 
integrating in v shows that 

div^; ^{A{v)ge) + 

= - V<1>(5,) 

+ fsiiZeG^^ + effe)) 

^-pV4>(x) + /, 

in the sense of distributions as e —)> 0, with the notation 


A{v) ■.= v^^-^\v\\ 


Using (2.12), we recast the limit above as 

divj, ^{A{v)ge) + Vg^) {9 - 0)V$ + fs 


(2.13) 


The elementary properties of the tensor field A used in the paper are re¬ 
called in the Appendix. In particular, elementary computations and Hilbert’s 
Lemma [T] show that 


A G (Ker£)-*- = Ran(£), 

so that there exists a unique tensor field denoted A such that 


CA = A and AT Ker C . 
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(See Lemma 0) Thus 

^{Mv)ge) = (^A{v)^Cge'^ . 

Returning to the scaled Boltzmann equation, we observe that 

=Q{ge,ge) - V ■ Vxge 
+ eM ^Va;<h • S/xiMge) 

- eM-^f,ix) • + eg,)) 

-^Q{g,g) - V ■ 

as e —>■ 0. 

At this point, we recall the following useful result. 

Lemma 2. For each </>,'!/’£ KerT, one has 

Q{4>,'ip) = \^^{H) ■ 


See formula (60) in [3] (one should notice the slightly different definitions 
of C and Q in formula (20) of |3], which account for the different sign and 
normalizing factor ^). 

Since g{x,-) G Ker£, Lemmaimplies that Q{g,g) = ^C{g‘^), so that 


On the other hand, (2.9) and (2.12) imply that 


g = 6 + u ■ V + 9^{\vf — 5) 


(2.14) 


so that 

g^ =A{u) ; A{v) + 26u ■ B{v) + + \0‘^{\v\^ — 5 )^ 

+ f + 2eu-v + ee{\v\^ -b), 

while 

V ■ = A{v) : VxU + B{v) ■ S/^O , 
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since u is divergence-free. Here, the notation B{v) designates the vector field 

B{v) := - 5)n, 

whose properties are recalled in the Appendix. 

Therefore 

^{{\A{v)g^ - A{v)v ■ V^t?)) 

=A(u) ; {\A{vr^) 

— Vu : {A{v) ® A{v)) 

=A{u) - u{Vu + {Vuf ), 

I and of the identities 

{AijA}^l) — ^ik^jl T ^^ij^kl j 

{AijA}^l') — k'l^SikSjl T Sil6jk ^^ij^kl) j 

with 

iy:=^(A:A}. (2.15) 

(See statement (2) in Lemma |^) 

Therefore 

div(A(n) - ia(Vu + (Vuf)) +Vq = (0- 0)V$ + /, , 
which can be recast as 


on account of (2.10 


div(n (g) u) — iaAu + Vp = 0V^(x) -j- fs{x) 


where p = q — This is precisely the motion equation in the 

Navier-Stokes-Fourier system. Notice that the term 


1 

e 




appearing on the left hand side of the equality (2.13) does not converge in 
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the sense of distributions in general, but its gradient does. Define 

T := lim V-(i|uPg'e). 

For each compactly supported, divergence free test vector held = £,{x), one 
has 

By Theorem 17’ in mi, T viewed as a 1-current is homologous to 0, which 
means precisely that T = Vq for some distribution q. 

Step 4- the heat conduction equation 

Next we explain how to derive the heat conduction equation. Mutiplying 
both sides of the scaled Boltzmann equation by 5) and integrating 

in V shows that 

diVa; ^{B{v)ge) = - Vd>(x) • (vg^) 

+ fs{x) ■ -I- ege)) 

—)• — u • V‘h(x). 

On the other hand 

-^{B{v) {\C{g'^) - V ■ VxQ)) 

={\B{v)g^ - B{v)v ■ Vxg) 

=(B(u)®2) . ud 

- {B{v)®B{v))-Ve 
=lu9-Kve. 
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Indeed, 

r , 

with 


k = {B-B). 


(2.16) 


(See statement (1) in Lemma |^) 

Therefore 

I div(u0) — kA0 = —u • Vd), 

which is precisely the heat conduction equation in the Navier-Stokes-Fourier 
system. □ 


2.2. Prom Kinetic Theory to Viscous Heating 

In the asymptotic Navier-Stokes-Fourier regime discussed above, the 
fluctuations of velocity field and of temperature are small and of the same 
order 0(e). In this case, the fluctuation of kinetic energy is negligible when 
compared to the fluctuation of internal energy. In order to keep both fluc¬ 
tuations small and of the same order, it is natural to scale the fluctuation 
of velocity field as 0(e), while the fluctuation of temperature should be of 
order O(e^). 

At the level of the Boltzmann equation, this scaling assumption is ob¬ 
tained by choosing the distribution function of the form 

Fe{x,v) = M{1 + eg^{x,v) + e^he{x,v)) (2-17) 


where 


ge{x,v) = —ge{x,—v), while he{x,v) = he{x,—v) (2-18) 

for a.e. {x,v). Here, we assume that there is no conservative force, i.e. we 
take the potential <I> identically 0. The total external force acting on the gas 
is therefore fs = /s(t) such that div/^ = 0. 
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The dimensionless Boltzmann equation (1.1) is scaled as follows: 


Ma = Ki = e, and Fr = 1. 


(2.19) 


In other words, the scaled Boltzmann equation takes the form 


V ■ + e^fs{x) ■ 


1 




( 2 . 20 ) 


e 


The exposition in this section follows closely [8], where the idea of the 
even-odd decomposition of the distribution function seems to have been used 
for the first time. 

First, we express the local balance laws of mass, momentum and energy 
in terms of the fluctuations and /ig. The odd contributions of either or 
vanish after integration in u, so that 



while 


V (g) vF^dv = / + e^(u (g) vhe) 



Hence, the local balance laws of mass, momentum and energy implied by 
the Boltzmann equation take the form 


div 3 ,(u 5 £) = 0 , (mass) 


while 


diva;(u (g) vhe) = fs + (^fs{he) , (momentum) 


and 


diYx{v\\v\^ge) = e^fs ■ {vg^) ■ (energy) 


Likewise, both sides of the Boltzmann equation are decomposed into 
even and odd components, observing that, for each rapidly decaying F 

C{Fo R) =C{F)o R, for all i? G 03 (R) • 
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Since the Maxwellian M is a radial function, one has M o R = M, and 
therefore 


C{(j) o R) = {Ccj}) o R , Q[cj) o R,if; o R) = Q{4>, o R 

for all R G 03 (R) and all rapidly decaying (j), ip. These identities are satisfied 
in particular for R = —I. Hence the even and odd components of C{F^) are 
respectively 

^(^^)even ^ ^ e^MQ{g,, g,) + e^MQ{K, h,) , 

C(T,)odd ^ -eMCg, + 2e^MQ{g,, K). 

Therefore, the scaled Boltzmann equation is equivalent to the system 

V • VxQe + e^M“Vs • '^viMge) = -Che + Q{ge, ge) + e^Q{he, he) , 

e\ • V^he - e^v • fs + e^M-^fs • V,(Mh,) = -Cge + 2e^Q{ge, he) . 

( 2 . 21 ) 

Theorem 2. Let Fe be a family of solutions of the scaled Boltzmann equation 
(2.20), whose relative fluctuations ge,he defined in {2.17)-(2.18) satisfy 

ge^ g, he^h 

weakly in x R^; (1 + \v\'^)Mdvdx), and 


Q{ 9 e,ge) Q{g, g), Q{ge, he) Q{g, h) 

weakly in L^(R^; (1 + \v\^)Mdv)), while 

e^iheVfiiv)) 0 

weakly in L‘^{VL) for each G L^(R^; (1 + \v\‘^)Mdv). 

Then 


g{x,v) = u{x) -v, 


( 2 . 22 ) 


while 


h{x,v) =\A{u{x)) : A{v) — A{v) : V xu{x) 

+ p{x) + {0{x) + i|n(x)p)l(|np - 3) 


(2.23) 
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where {u, 9) is a solution of the incompressible Navier-Stokes-Fourier system 
with viscous heating 

div u = 0 , 


div(M®^) + Vp = v/S.u + fs , 


(2.24) 


\T|2 


I div(u0) — u ■ VxP = + ^ly \ VxU + (Vxu)^ 

The values of the viscosity v and heat diffusivity k are determined implic¬ 


itly in terms of the collision integral, by formulas (2.15) and (2.16), as in 
Theorem [ij 


Proof. The proof follows more or less the same lines as that of Theorem 
see also [8]. 

Step 1: asymptotic form of g^ and divergence-free condition 


We first deduce from the second equation in (2.21) and the assumption 
that —)• 5 in the sense of distributions that 

£5 = 0. 

Hence g{x, •) G Ker £ and v i-A g{x, •) is odd for a.e. x, so that g is of the 
form (2.22[). The local conservation of mass implies that 


0 = Ai\x{vge) div3;(ug') = diVxU 


in the sense of distributions, so that 


divi; u = 0 . 


(2.25) 


Step 2: asymptotic form of h^ and divergence-free condition 


Next we deduce from the first equation in (2.21) that 


he ^ h with V ■ Vxd = —£/i + Qig, g) ■ 
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Since g{x,-) G Ker £, applying Lemmashows that 


C{h — 29 ) — ■ ^x9 = —V ® V : = —^iv) ■ '^xU, 


since u is divergence free by (2.25). Hence 


h{x, v) = y {x, v) - A{v) : VxU + /ro(x, v ), 
with /io(T, •) G Ker/1 for a.e. x. Since /iq is even in v, it is of the form 


/io(x, v) = Tu{x) + 0(x) 2 (|nf — 3). (2.26) 


On the other hand 

g‘^{x, v) = {u{x) ■ v)"^ = u{x) ® u{x) : v 
= A{u{x)) : A{v) + \\u{x)\^\v\^ , 

so that h is of the form 

h{x, v) =lA{uix)) : A{v) - A{v) : Vxu(x) 

+ Pix) + {0{x) + 5 |n(x)p)|(|n|^ - 3), 

with p = w + g|np. 

Step 3: motion equation 

Passing to the limit in the sense of distributions in the local balance law 
of momentum shows that 


diva:(n (g) vh) = fs . 

We insert the expression for h found above in this identity. Observe that 
{v ® vh) =(H(n)(^H(u(x)) : A{v) — A{v) : Vxu(x))} 

+ (P(^) + (^(^) + lHx)\‘^) - 3 )))-^ 

because A and A _L Ker/1. By statements (l)-(2) in Lemma 
{A{v){lA{u{x)) : A{v) - A{v) : V^u(x))} = A(u(x)) - u{VxU + (Vxu )'^), 
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while 


(llwp {p{x) + (6»(x) + ^\u{x)\‘^)l{\v\‘^ - 3))) = p{x) + e{x) + 5 |ii(x)p 


Hence 


diva: {A{u) - v{VxU + {Vxuf)) + V^{p + 9+ = fs , 


or equivalently 


diVa;(ii ®u) - vA^u + VxP = fs 


(2.28) 


with p = p + 9. 


Step 4-1 heat equation 

Finally, we combine the local balance of mass and energy so that 

diva, \ {Bge) = fs ■ {vge) fs-u 

in the sense of distributions. 

Next, we transform the left hand side of the equality above by exactly 
the same method as in the previous section; 

^{Bgf) = l^B^Cg^j 

= {mQ{g„ K)-v - e^M-^fs • V„(Mh,))) 

^ {B{2Q{g,h)-vV,h)). 


Notice that 


/R3 


because B _L Ker£. Hence 


Bfs ■ VyMdv = -{Bfs ■v) = 0, 


dwx{B{2Q{g, h) - v ■ Vxh)) = fs ■ u , 


and we insert in this last expression the explicit formulas for g and h. 
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First 


22(5, h) =2Q{u ■ V, lA{u) : A{v) + p + {9 + l\u\‘^)l{\v\'^ - 3)) 
— 2Q(tt • V, A{v) : Vxu). 


(2.29) 


Applying formula (2.18c) in [8] shows that 

2Q{u ■ V, IA{u) : A{v) + p + {9 + ^\u\‘^)^{\v\‘^ - 3)) 

= C{9u ■ B{v) + ^C{v) : u), 

where 

C{v) := 0 V ® V — 3v ® I). 

Hence 

2{B{v)Q{u ■ V, lA{u) : A{v) + p + {9 + g|'u|^) 5 (|'Cp - 3))) 

= {B{v)C{9u ■ B{v) + lC{v) :u^u^u)) 

= {B{v){9u ■ B{v) + ^C{v) : tt (g) u (g) tt)) 

= 5 (|n|^ + 59)u 

by Lemma and statement (1) in Lemma 

By Proposition 2.6 in [H] and statements (1) and (3) of Lemma 

2{B{v)Q{u ■ V, A{v) : V^u)) ={{A{v) ■ u)A{v) : V^u) 

— {B{v){u ■ v)A{v) : Vxu) (2-31) 

={u - Ik){VxU +{Vxu)'^) ■ u . 


Finally, we compute 

{B{v) (g) vh) =^{B{v) (g) u (g) A{v)) : A{u) 

— {B{v) (g) u ® A{v)) : VxU 
+ {9 + l\u\‘^){B{v) (g) B{v)) 

= |kA(u) — c{VxU + (S/xu)'^) + k{9 + g|ttp)/ 

according to statements (1) and (3) of Lemma Notice that 

p{B (g) u) = 0 , 
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and that 


because B _L Ker£. 




Therefore, putting together (2.30)-(2.31)-(2.32), we arrive at the identity 


diVx{B{2Q{g, h) - v ■ V^h)) 

= diVa;( 5 (|up + 59)u - {v - ‘Ik){VxU + (V^u)^) • u) (2.33) 

-V,r (g) Va, : (|kA(u) - c{VxU + {Vxuf) + k{6 + ^\u\‘^)I). 

This identity can be substantially simplified, as follows. First, 

Vx ® Va; : (Va;U + (S/xu)'^) = 0 


because of the divergence-free condition (2.25). On the other hand. 


diVa; ((Va:U (VxU)'^) ■ u)) = Vx <Si Vx '■ U (gi U + Ax ^ 




while 

Va; ® Va: : (A(u) -I- ||up/) = Va; ® Va; : U ® U — ^Aa;|u|^ -|- |Aa;|u|^ 

= Vx <8> Va; : U U + Aa:2|'u|^ . 


Hence 


diVa;(.B(2Q(g', h) - V ■ Vxh)) 

= diva; (jdup -I- 56)u - V {VxU -|- (Vxu)'^) • u) - kAxO , 
so that the limiting form of the local energy balance is 


diVa: ( 2 ( 1 ^ 1 ^ + 56)u) = l^diVx{{VxU + (VxU)'^) ■ u) -I- kAx 9 + fs ■ u . 


On the other hand, multiplying both sides of the Navier-Stokes motion equa¬ 
tion by u, we arrive at the identity 

diVa: {u\\u\^^ + u ■ VxP = fs ' U + UU ■ AxU 

= fs-U + VU- diVa; (Va;U (Va^u)'^) . 
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We further simplify the right hand side of the equality above as follows: 


u ■ diVa; [VxU + (VxU)'^) = diVa; {{VxU + (Vxu)'^) u) 

- (Va:M + {VxU)'^) ■ VxU 
= diVa; ( (Va;U + (Va;U)^) u) 

- h |Va:U + (Va;U)^|^ . 


Combining these identities with (2.33) above leads to 


I diVa;(u0) — U ■ VxP = kAx 8 + ll' |Va;M + (Va;tt)'^|^ . 


(2.34) 

□ 


Observe that the Navier-Stokes motion equation is exactly the same 


when <1> = 0 for both scaling assumptions (2.1) and (2.19). The temperature 


equation, however, is very different according to whether the Froude number 


is 0{^/e) as in (2.1), or 0(1) as in (2.19). 

The viscous heating term on the right hand side of the equation gov¬ 
erning the temperature field appears in Sone’s asymptotic analysis of the 
hydrodynamic limits of the Boltzmann equation in the weakly nonlinear 
regime — see [l6] . Sone’s original work [Uj on the weakly nonlinear hydro- 
dynamic limits of kinetic theory was written in the case of the BGK model; 
see IE! for the extension to the Boltzmann equation. Sone’s argument is 
based on the Hilbert expansion. One should pay attention to a particu¬ 


lar feature of Sone’s theory: the pressure field p in (2.8) and the velocity 


and temperature fields u and 9 in (2.8) do not appear at the same order in 
Sone’s expansion; in fact p appears at order O(e^) in the expansion of the 
distribution function in powers of e, while u and 9 appear at order 0(e) in 
that same expansion. See formulas (3.77), (3.79b-c), (3.80d) and (3.88a-c) 
in section 3.2.2 of [E]. The limit leading to (2.24) corresponds to a situation 


where u appears at order 0(e) in Sone’s expansion, while the leading order 
temperature fluctuation (denoted tsi in Sone’s analysis, see formula (3.79c) 
in m) is identically zero. The temperature fluctuation appears at order 
O(e^) in Sone’s expansion, together with the pressure field p, and the tem¬ 


perature equation in (2.24) coincides with formula (3.89c) in section 3.2.2 
of [E]. (Sone’s analysis in |E] does not involve an external force, but the 
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work of the external force disappears from the temperature equation when 
combining the motion equation and the energy equation as explained above.) 


2.3. Boundary Conditions 


The discussion of the incompressible Navier-Stokes-Fourier limit of the 
Boltzmann equation presented above would remain incomplete without dis¬ 
cussing the boundary condition. In this section, we briefly describe the 
simplest imaginable situation. 


Assume that the scaled Boltzmann equation (2.2) is supplemented with a 
diffuse reflection condition at the boundary of the spatial domain 12 on which 
the Boltzmann equation (1.1) is posed. In other words, for each x G 512, one 
has 


F^{x,v) = V^M / Ff:{x,v){v ■ nx)-\-dv , v ■ Ux < 0 , 
JR3 


(2.35) 


where x i—)• is the unit normal field defined on the boundary 512 of the 
spatial domain. Here, we have assumed for simplicity that there is no tem¬ 
perature gradient on 512. The constant temperature at the boundary (i.e. 
the temprature 1 in the Maxwellian state M) defines the scale of the speed 
of sound in the interior of the domain. 

By construction 


/R3 


F^{x, v)v ■ Uxdv = 0 , X G 512 , 


which means that the net mass flux at each point x G 512 is identically 0. 


This suggests that the boundary condition (2.35) should be supplemented 
with the additional condition 


QxR3 


F^{x, v)dxdv = |12| , 


(2.36) 


that is consistent to leading order with the normalization of in (2.3) 


(2.17), and is equivalent to the condition (2.4) already introduced above in 
the case Ma = Ki = Fr^ = e, and 


{he)dx = 0 
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in the case Aia = Ki = Fr = e. (Notice that, in the latter case, (g^) = 0 a.e. 
on 0 since is odd in ?;.) 

Besides, we assume that the force fg satisfies both 

divj; fs = 0 on n, fs'nx = 0 on dO,. (2.37) 


Define 

Axicp) := V^{4>iv ■ nx)+) ■ 


Theorem 3. Let Ma = Kn = Fr^ = e, and consider a family of solutions of 
the scaled Boltzmann equation (2.2) supplemented with the dijfuse reflection 
condition (2.35) and with the total mass condition (2.36). Assume that 
Fg = M{Zge'^^ + egflj as in (2.3) and satisfies the same assumptions as in 
Theorem(^ Assume moreover that the family of traces of g^ on the boundary 
dLl X satisfies 


S'laoxRS weakly in L^{dVt x R ^;\v ■ nx\MdvdS{x)) 
where g is such that g^ ^ g weakly in L^(D x Mdvdx). Then 
g{x, v) = 0 + u{x) ■ V + 9{x)\{\v\^ — 5) 


where 




and {u,9) is a solution of the Navier-Stokes-Fourier system (2.8), with Diri- 
chlet boundary condition 


u\ 


an 




Proof. Indeed, the diffuse reflection condition implies that 
gflx,v) = Axigeix,-)), v-nx<0. 


Thus one can pass to the limit as e —>■ 0 in (2.38). One arrives at 
g{x, v) = Ax{g{x, •)), x G 5D , u G R^ . 


(2.38) 
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Since we already know from Theorem that g of the form 
g{x,v) = 9 + u{x) ■ V + 6{x)\{\v\^ — b) ^ 


this implies that 




= 0. 


□ 


Theorem 4. Let Ma = iih = e while Fr=l, and consider a family of solu¬ 


tions of the scaled Boltzmann equation (2.20) supplemented with the diffuse 


reflection condition (2.35) and with the total mass condition (2.36). Assume 


that Ff: satisfies the same assumptions as in Theorem and that the odd 
and even part of the relative fluctuation of distribution function, resp. g^ 


and hf. defined in (2.1'}) are continuous in v and satisfy the condition 




^G9QxR3 


l(9nxR3 


locally uniformly in x,v, where we recall that g and h are the weak limits of 
g^ and in x R^; (1 + \v\‘^)Mdvdx) as e —>■ 0. Then g and h are given 
by the expressions (2.22) and 12.25), where {u,9) is a solution of the system 


(2.24) with the Dirichlet boundary condition 

^lan ^ ® • 


^l(9Q “ ® ■ 


Proof. Specializing the equality above to the case where v is tangential to 
the boundary, we find that 

gflx, v) = K{{ge + ehf){x, •)) - ehflx, v), x G dQ , v ■ = 0 , 

and observe that the left hand side of this equality is odd in u = v — {v-nx)nx, 
while the right hand side is even. Therefore both sides vanish, so that 

gflx, u) = 0 , X G dPl , V ■ Ux = 0. 


while 


he{x,v) 


Ax 



9e + h, 


(x,-) 



X G dLl , V ■ Ux = 0 . 
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Passing to the limit on both sides of the first equality as e —)• 0 shows that 


g{x, v) = Kx{g{x, ■)), X G dQ, V ■nx = 0, 


and we conclude from (2.22) that 




Next we consider the differential operator 

P(x, Dv) •— ^ V 

— which is the orthogonal projection of V„ on the tangential direction of 
dn at X, and observe that 

P{x, Dx)hf:{x, v) = 0 , X G dQ , v ■ Ux = 0 . 

Passing to the limit in both sides of this identity, we conclude that 

P{x, Dx)h{x, v) = 0 , X G dCl , v ■ Ux = 0 ■ 


Substituting the expression (2.27) in this equality, we find that 


P{x, Dy)h{x, v) = —P{x, Dy)A{v) : Vu{x) + Vt9{x) , x G di} , v ■ Ux = 0 


with 


Observe that 


Vr = v - {v ■ nx)nx 


P{x, Dy)A{v) : Vu{x) = Vr ■ Vtt(x) + {I — Ux 0 nx)V{u{x) ■ Vr) = 0 , 

since we already know that = 0 and all the derivatives of u appearing 
in the expression above are taken in directions tangential to di}. Hence 

Vr9{x) = 0 , for each x G dCl , v ■ Ux = 0 , 


which implies that 9\g^ = 0. 


□ 
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3. Spatially Periodic Steady Solutions 
3.1. The Navier-Stokes Equations 

In this section, we assume that the spatial domain is = T^. Consider 
the system (2.8) posed on T^, and seek solutions satisfying 


/ u{x)dx = / 0{x)dx = O. 

Jt3 

For simplicity, we assume further that = 0. 


(3.1) 


Multiplying both sides of the last equation in (2.8) by 9, we see that 


u ■ V(l 02 ) = div(ul 02 ) = K 0 A 0 , 


so that 

K [ \V9{x)\‘^dx = 0 . 

JT3 

Therefore 0 = 0. 

Conversely, if n is a solution of the motion equation and ‘h = 0, then 
{u, 0) is a solution of the Navier-Stokes-Fourier system without viscous heat¬ 
ing. 

There exist indeed nontrivial solutions of the motion equation with 
nonzero solenoidal external force fg. The simplest example is the case of 
a shear flow 


u{x) = (0,0, f7(xi,X2)), fs{x) = (0,0, a(xi,X 2 )) 


with 


Obviously 



U{xi, X 2 )dxidx 2 



a{xi,X2)dxidx2 = 0 . 


div u = dx^U = 0 , div fg = dx^a = 0 , 


and 


div(u®^) = {u ■ V)n = Udx^u = 0 , 
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SO that the Navier-Stokes equation reduces to the Poisson equation in T^: 
-i'AU{xi,X2) = a{xi,X2) , ixi,X2) G . 


For each zero-mean a G L^(T^), there exists a unique zero-mean U G 
of the Poisson equation above. 

More generally, the following result is classical. 


Theorem 5. For each fg G L^(T^) satisfying 



fs{x)dx = 0, 


there exists at least one solution u G of the Navier-Stokes equations 

with external force fs such that 



u{x)dx = 0. 


Besides, there exists c > 0 such that the solution u G and is unique 

if ||/s||l2 < CZ22. 


The proof of this classical result is given below — see chapter II, §1 in 
|49| for a similar result in a slightly different (nonperiodic) setting. First, we 
recall some elements of notation. We denote by the subspace of L^(T^; R^) 
of vector fields v such that 

/ v{x)dx = 0 , 

7t3 

and we set = R^(T^) n S^. It will be convenient to use the norm 

||u||^i := \\Vv\\l2 . 

We denote by Ft the L^-orthogonal projection on divergence-free vector 
fields. In other words, if u G ij, its Fourier decomposition is 

V{X)= Y. 

fcez3\{o} 
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and 

Iiv{x)= ^ (ufc - . 

fcGZ3\{0} 

Likewise, for each zero-mean (p G Sj we define 

(—A)“®(/)(x) ;= where (f>k = f dx . 

fcez3\{o} 

Proof. Consider the map 

T: (-A)-indiv(r;®2). 

First, observe that T maps into itself: 

\\Tv\\Hr = ||(-A)-V2ndiv(a®2)||^, = ||n(-A)-V2div(z;^2)b. 

< ||(-A)-V2div(z;®2)||^^ < ||^®2||^^ < ||^||2^ < ^2||^||2^ 

by Sobolev embedding C L^(T^) for 1 < p < 6). Similarly 

\\Tu - Tv\\hi = ||(-A)-^/2ndiv(M®2 _ v^‘^)\\l2 < - V^^\\l2 

< ||ii (g) (tt - f)||2,2 -I- II (ti - v) (g) tj||x,2 < (||rt||i,4 -I- 11 ^ 111 , 4 ) ||ii - tj||2,4 

^ c'^dl^lUi + IklUOII'*^ “ • 

Hence T is continuous from into itself, and Lipschitz continuous on balls 
of 

The Navier-Stokes equations can be put in the form 
u = {-uA)-^fs - iTu, 
and is embedded into the family of equations 

Ta{u) = {-uA)~^fs , where Ta = I + 
parametrized by a G [0,1]. 

If u, u G with ||u||/^i < R and ||u||j|^i < R, then the map Qa defined 
by 

Ga(w) = {-uA)~^fs + W- Rol{w) 
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satisfies the bound 

2C^ R 

WGaiu) - Gaiv)\\H^ < ^ ||u - v\\hi 

for a G [0,1]. On the other hand, if u G with ||u||/^i < R 
\\Qa{u)\\H^ < ||(-l/2A)-V2/,b2 + <R. 

Thus, if 

||(_A)-V2j^||^ 2 < andi?<^, 

then Qa maps the closed ball B{0,R) C S)i into itself and is a strict con¬ 
traction on 5(0, i?). Hence Qa has a unique fixed point in 5(0,5) C for 
a G [0,1] provided that 

||(_A)-V2j^||^ 2 < ^ , with R = ^||(-A)“^/Vs|U 2 • 

In particular, for a = 1, this unique fixed point of Qa is the unique solution 
of the Navier-Stokes equation in 5(0, R) C 

The estimate 

\\Tu - Tv\\h^ < {\\u\\l4 + ||u||£,4)||u - v\\l4 

and Rellich’s theorem imply that the map T is compact in Indeed, 
if ttn —?• u weakly in then Un ^ u strongly in L^(T^) by the Rellich 
compactness theorem, and the inequality above with v = Un shows that 
Tun ^ Tu strongly in On the other hand, for each u £ the equation 
Tq(u) = {—vA)~^fs is equivalent to 

div u = 0 , a{u ■ V)u -|- Vp = vAu + fs, x G T^ . 

Multiplying both sides of the motion equation by u and integrating over T^, 
one finds that 
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SO that 


\\Vuh2<l\\{-A)-y^fsh2. 


Setting R' = ^||(—A)“^/^/s||i 2 + 1, we see that Ta maps B{0,R') C into 
and that u G dB{0,R') C S)^ implies that Fa{y) 7 ^ {—vA)~^fs for all 
a G [0,1]. Therefore 


degiee{Ri,B{0,R'),{-iyA) ^ f^) = degiee{Ro, B{0, R'), {-i^A) V^) = 1 
because Rq is the identity. Hence the equation 

^i{u) = (-z^A)“Vs, 

which is equivalent to the steady Navier-Stokes equations in T^, has at least 
one solution in □ 

However, the only solution of the Navier-Stokes-Fourier system with 


viscous heating (2.24) satisfying (3.1) is (0,0). Indeed, integrating both 


sides of the last equation in (2.24) shows that 


/ |Vri + (Vn)^|^(x)dx = 0 . 
Jt3 


Hence 


Vu + (Vu)^ = 0 . 

In particular, div u = 0 and 

— div(VM) = —Au = div((VM)'^) = V(divM) = 0 , 


so that tt is a harmonic vector field on T'* satisfying (3.1). Hence u = 0 


Returning to the heat equation in (2.24), we see that 0 is a harmonic function 
on T^, and (3.1) implies that 9 = 0. 
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3.2. The Boltzmann Equation 

Theorem 6. Let F = F{x,v) be a solution of the steady Boltzmann 
equation 


V ■ VxF + f{x) ■ VvF = C{F), X ^ , V ^ , 

and assume that F is rapidly decaying in v while In F has polynomial growth 
in V as |z;| —)■ oo. 

Then f is a gradient field, while F is a Maxwellian distribution with 
constant temperature. More precisely, there exists G C'^(T^), a vector 
u G and two constants 6 > 0 and C > 0 such that 

f = — V<1> , and F{x, v) = C exp (sl'^ “ 

with 

u ■ V$ = 0. 


In particular, if div f = 0, then / = 0 and F is a uniform Maxwellian. 


Proof. The global form of Boltzmann’s H theorem shows that 


C{F) In Fdxdv = 0 . 

/T3xR3 

Hence F(x,v) is a local Maxwellian satisfying 

V ■ VxF + f{x) ■ VvF = C{F) = 0 , x £ , V £ R^ . 

Setting F{x,v) = M.(p(x),u(x),e{x))iv), the Boltzmann equation reduces to 
v-Vp{x) (v — u) 0 V : Vu{x) v-V6{x) 


p{x) 


+ 


9{x) 


+ 


2e{x) 


-(|u —u(x)| —39{x)) 


= fix) 


V — u{x) 

9{x) 


Setting 


V := 


V — u{x) 

VW) 
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the equality above is recast as 

u{x) • V in p{x) + 6{x)~^^^V • (0(x)V in p{x) + {u • V)u) 

m\Vf - 3)u(x) • Vln0(x) + 

+(|yp - 3)y • v^/^ = e{x)-^/^f{x) ■ v. 

Hence 

vV0 = 0 


so that 


0 = Const. 


and the equality above reduces to 

u(x) -Vlnpix) + . (0vin/)(x) + (u ■ V)u) + 

= e-^/^f{x)-v. 

Therefore 

u ■ Vp = 0 , 
eVlnp+{u-V)u = f, 

Vu + (Vu)'^ = 0 . 

Using the third equation, the second equation is recast as 

eV In p+{u- V)u = eV\np- {Vufu = 9V\np- V\\u\^ = f 
so that 

/ = V{einp- l\u\^) 

must be a gradient field. 

Next 

Vu + (Vu)"^ = 0 

which implies as above that u is harmonic on T^, and therefore is a constant. 
We conclude that 

/(x) = —V<f>(x), F(x, u) = with u • V<h(x) = 0 . 
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Hence u = 0 unless takes its values in an affine space orthogonal to u. □ 


4. Physical discussion and study of a numerical example 

As we observed in Sec. 3.2, there is no spatially periodic steady solution 
of the Boltzmann equation for an external force that is not derived from a 
potential. In contrast, the Navier-Stokes equations have a spatially periodic 
and steady solution for such an external force. This discrepancy seems to 
be contradictory, since the Navier-Stokes equations are derived from the 
Boltzmann equation, as shown in Sec. 2.1. In this section, we will further 
examine this seemingly contradictory results. 


4.1 A possible physical explanation of the paradox 

Let us consider a gas in a periodic box or in a specularly reflecting box. 
As we have seen in Sec. 3.2, the steady solution of the Boltzmann equation 
with an external force derived from a potential has the following properties: 
the temperature is uniform, the macroscopic flow of the gas vanishes except 
for a special case, and the density is distributed according to the potential 
(i.e., a stratified gas at rest). Now, let us consider a time-dependent problem 
starting from a given initial state. If the external force has a potential, 
the time-dependent solution should approach the steady solution mentioned 
above, i.e., the solution without a gas flow and with a density stratification, 
in the long-time limit. Then, what will happen when the force does not have 
a potential? At the initial stage, a gas flow is caused by the external force. 
But, since the force does not have a potential, the density stratification that 
blocks the gas flow cannot be formed. Therefore, the flow remains forever, or 
at least, for much longer time. The induced gas flow, in general, has a shear. 
If we consider the case with relatively small Knudsen number, this shear 
gives rise to the viscous heating. However, because the boundary is periodic 
or adiabatic (in the case of specularly reflecting box), the heat generated by 
the shear in the gas cannot escape through the boundary. This means that 
the temperature in the gas will increase indefinitely. This is the reason why 
there is no steady solution for the Boltzmann equation. 
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On the other hand, the Navier-Stokes equations (2.8) are derived in the 
limit where the effect of the viscous heating is negligibly small. Therefore, 
the mass and momentum equations (the first two equations in (2.8)) are 
decoupled with the energy equation (the last equation) when $ = 0. The 
former equations give a velocity field with a shear, and the latter equation 
gives a constant temperature field. However, if we include the effect of 
viscous heating, as was done in Sec. 2.2, the energy equation is changed 
to the last equation in (2.24), the mass and momentum equations being 
unchanged. Therefore, the flow velocity field is the same, but the new energy 
equation does not have a solution for this velocity field. If we consider a 
time-dependent version of this energy equation, we can easily see that the 
temperature increases indefinitely. In conclusion, the difference between the 
Boltzmann equation and the Navier-Stokes equations is due to the fact that 
the effect of viscous heating is neglected in the Navier-Stokes limit with 
Ma = Ki = Fr^ = e. 


4.2 Numerical example 

For the purpose of understanding the phenomenon of heating predicted 
above, we consider a simple numerical example. Let us consider a gas in 
a two-dimensional square box —1/2 < xi < 1/2, —1/2 < X 2 < 1/2 with 
periodic condition on each side. We assume that the external force is of the 
form / = (0, /osin27rxi, 0), which is divergence free and does not have a 
potential. Initially, the gas is in a uniform equilibrium state at rest with 
density 1 and temperature 1. We pursue the time evolution of the solution 
and observe whether the temperature increases indehnitely or not. 

We analyze the problem mainly using the Bhatnagar-Gross-Krook (also 
known as BGK) model, but some results based on the original Boltzmann 
equation will also be presented. Because of the form of the external force, 
we can assume that the flow field is spatially one dimensional depending 
only on xi and is periodic in xi with period 1. Moreover, the external 
force is symmetric with respect to xi = —1/4 and xi = 1/4, so that we 
can also assume the same symmetry for the flow field. Therefore, placing 
specularly reflecting boundaries at xi = —1/4 and = 1/4, we can analyze 
the problem in the hnite interval —1/4 < xi < 1/4. 
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Here, we formulate the problem using the BGK model. In the present 
problem, the BGK model in an appropriate dimensionless form can be writ¬ 
ten as 

dtF + vid^^F + /o sui2TTXid^^F = 

where K(xi, t, v) is the velocity distribution function including the time vari¬ 
able t, and 


p(M-F), 


Cbgk(F) = 

p(xi,t)= / Fdv, 

Ui(xi,t} = ViFdv, 

pJn^ 

0{xi,t) = ^ I \v — u\‘^Fdv. 


pixi,t) 


[27r0(xi, 


exp 


(i = 1, 2, 3; U3 = 0) 


V - u(xi,t)p \ 

26'(xi,t) ) ’ 


Here, Kn is the Knudsen number, i.e., the mean free path in the initial 
equilibrium state at rest divided by the length of the period (note that 
Kn is denoted by e in Sec. 2 because the limit as e —>• 0 is discussed there). 
The factor (S/vr)^/^ appears because of the manner of nondimensionalization 
used here in consistency with the form of the Boltzmann equation in earlier 
sections. The specular reflection condition at xi = —1/4 and 1/4 is as 
follows: 


F(±l/4, t, x) = F(±l/4, t, for =p xi > 0, 

where R is the reflection operator: Rv = (—xi,X 2 ,X 3 ). The initial condition 
is given by 


F(xi,0,x) = A4(i^o,i)- 

We solve this initial-boundary-value problem by the finite-difference method. 

Now we show some of the numerical results for Kn = 0.1 and for /o = 2 
(Fig. 4.1) and 0.2 (Fig. 4.2). 
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Figure [4yT| a) shows the profile of the X 2 component of the flow velocity 
U 2 from f = 0 to 1000/\/2 in the half interval 0 < < 1/4. Since ui and U 2 

are odd functions of xi , and p and 9 are even functions of xi , we show the pro- 
hles of these quantities only in the half interval here and in what follows. The 
sinusoidal external force parallel to the X 2 direction induces U 2 at the very 
early stage, but U 2 starts decreasing after t = 1.5/\/2. Figure 4.1'b) shows 


the time evolution of the average temperature 9av{t) = 2 f^^-^^^9(xi,t)dxi 

c) the corresponding evolution of the profile 


4.1 


until t = 1000/\/2, and Fig. 
of the deviation 0 — 9av The 0 becomes nonuniform at the early stage but 
tends to become uniform after t = 1.5/\/2. On the other hand, Oav increases 
and reaches 50 at t = 1000/-v/2- As shown in Fig. the density p is 

nonuniform only at the very early stage and is almost uniform (i.e., almost 
p = 1) at t = 30/-v/2. Corresponding to the nonuniformity of the density, 
the flow-velocity component ui perpendicular to the external force arises at 
the very early stage [Fig. |4.1[e)], but its magnitude is very small and prac¬ 


tically vanishes at t = lt)/^/2. Figures 4.1 f) and |4.l| [g) show the long-time 
behavior of 6av{t) and the average speed, |u 2 |ai)(i) = 2 |u 2 (xi, t)|(ixi, 

up to t = 24000/-v/2. The left hgures show the double-logarithmic plot of 6av 
versus t and that of |M 2 |a'!; versus t. In the right figures, the gradients of the 
curves in the left hgures, i.e., a = dln9av/dlnt and /3 = d\Ti\u 2 \av/dtat are 
plotted. If a and /3 approach constant values, say ao and /3o, respectively, 
then we have the long-time behavior as Oav ~ and |rt 2 |a'u ~ Cut^° with 

positive constants Cg and Cu- From Figs. |4.1| [f) and 4.1 [g), it is still not 
clear whether a and /3 converge to hnite values or not. But, if it is the case, 
it is likely that ao ~ 0.66 and /Iq ~ —0.34 ps — (1 — ao). 

Figures [T^a)-(T2|[g) show the behavior, corresponding to Figs. |4.l)[a)- 


4.1 [g), for a weaker external force (/o = 0.2). The tendency of the time 


evolution of the solution is similar to Fig. |4.1[ However, since the magnitude 
of the force is 1/10, the resulting flow and the temperature rise are smaller. 
As Fig. |4.2| [a) shows, the how speed |tt 2 |, which is smaller by one order of 
magnitude, takes the maximum at around t = ^1^/2 and decreases more 
slowly than in Fig. |4.1| [a). Figure |4.2| [b) shows that the increase of \9av\ 
is much slower compared to Fig. |4.1| [b). One sees from Figs. |4.2[ c)- |43| [e) 
that the nonuniformity of 9 and p and the magnitude of ui are smaller by 
two orders of magnitude. In Figs. 4.2 f) and 4.2 [g), we show the long-time 
behavior of 9av and |tt 2 | up to an extremely large time, t = 336000/-v/2. As 
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in 


Figs. E3f) and |4.1[ g), the left figures are the log —log plots, and the 
right figures are their gradients a and j3. Even at such a large time, it is not 
clear whether or not a and /3 converge to constants. But, if they converge, 
the values would not differ much from the case of /o = 2 [cf. Figs. |4.1[f) and 


Finally, we present some results based on the Boltzmann equation. Let 
us replace the right-hand side, (l/Kn)CBGK(.F), of the basic equation with 
(l/Kn)C(E), where C is the dimensionless Boltzmann collision operator for 
hard-sphere molecules, as in Sec. 1. The mean free path used here to define 
Ki is (\/27r(i^no)~^ with d the diameter of a molecule and no the refer¬ 
ence molecular number density, which is related to £o in Sec. 1 [cf. 0] as 
(2\/27r)“^^o- We employ the direct simulation Monte Carlo (DSMC) method 
as the solution method. The method is a particle and stochastic one, so that 
it contains the inherent statistical fluctuations. For steady problems, we 
can take the time average over a long interval of time to reduce the fluctu¬ 
ations. But the time-averaging does not work for time-dependent problems 
as the present one. The only possible way to reduce them is to perform 
many independent runs and take an ensemble average over the runs. The 
method is also not appropriate for describing small quantities because they 
are hidden in the fluctuations. In fact, it is impossible to obtain 9 — 6av, p, 
and ui- Nevertheless, we show in Figs. 4.3 and 4.4 the results for the same 


cases as Figs. 4.1 and |4.2[ respectively. That is. Fig. 4.3 is for Kn = 0.1, 
/o = 2, and Fig. 4.4 for Kn = 0.1, /o = 0.2. Figures 4^a) and |4.4|[a) 


correspond respectively to Figs. |4.1[ a) and |4.2[ a); Figs. 4.3 [b) and |4.4[ b) to 
Figs. and |4.2| [b); Fig. |4.3| [c) to Fig. |4.1| [c) (it is impossible to obtain 


the corresponding figure for /o = 0.2); Figs. |4.3[ d) and |4.4[ c) to Figs. |4.l| (f) 
and |4.2[ f); and Figs. |4.3| (e) and |4.4[ d) to Figs. |4.l| (g) and |4.2K g). Figure 
4.3 shows the result of the ensemble average over 96 independent runs. In 


addition, in Fig. |4.3[ e), (3, which is the average of f3 over the time interval 
[t - 500/\/2, t], is shown instead of (3 itself, since it is impossible to obtain 
the plot of /3 in a reasonable form. In Fig. 4.4, we show the result based on 
the ensemble average over 96 independent runs for 0 < t < 100/\/2, that 
based on the ensemble average over 24 independent runs and the time aver¬ 
age over the interval [t — 5/\/2, t -|- 5/-v/2] for 100/-v/2 < t < 5000/\/2, and 
that based on the ensemble average over 12 independent runs and the time 
average over the interval [t — 10/\/2, t + 10/\/2] for 5000/\/2 < t. The /S in 
Fig. 4.4 d) is the average of /3 over the time interval [t — 10000/-v/2j t]. The 
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(a) U 2 (\/2i < 60000) 

, 15,-™- 



(c) e - (\/2t < 12000) 


Figure 4.2: BGK model with Kn = 0.1 and /o = 0.2 
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(a) U 2 (\/2i < 1000) 






(C) 0 - 0av (^21 < 30) 


Figure 4.3: Boltzmann equation with Kn = 0.1 and /o = 2. The ensemble average 
over 96 independent runs is shown. The /3{t) is the time average of (3 over the 
interval [t — 500/-\/2, t]. 
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0.35: .. . 
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°'^^00 800 1000 
{y/2t < 1000) 


Figure 4.3: Continued 


computed time in Fig. 4.3 (or Fig. 4.4) is much shorter than that in Fig. 4.1 


(or Fig. 4.2) because the longer computation is practically impossible. How¬ 
ever, the convergence of a seems to be faster for hard-sphere molecules. It 
is likely that a tends to approach about 0.665 in Fig. 4.3 )d) and about 0.66 
in Fig. |4.4[ c). Therefore, irrespective of the model of the collision term, we 
have the asymptotic behavior like Qav ~ for large t. 


4.3 Interpretation of the nnmerical results 

Next, in order to understand the slow increase of the temperature and 
the slow decrease of the flow speed, we try a rough discussion on the basis 
of the compressible Navier-Stokes equations. The numerical results show 
that the flow is almost unidirectional, i.e., u k. (0,U2,0). Therefore, let us 
assume that the flow is unidirectional, i.e., u = ( 0 , U 2 , 0 ) and the problem 
in spatially one-dimensional {d/dx 2 = d/dxs = 0). Then, the compressible 
Navier-Stokes equations reduce to the following equations: 
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(c) 0av (\/2t < 21000) 




15000 


20000 


Figure 4.4: Boltzmann equation with Kn = 0.1 and /o = 0.2. The result shown 
here is based on the ensemble average over N independent runs and time average 
over the interval [t — C 1^/2, t + 01^/2]^ where N and C are chosen appropriately 
depending on t. The details are given in the main text. The /? is the average of f3 
over the time interval [t — 10000/-\/2, t]. 
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^=0 

dxi 


du2 


1 d 

p dxi 



+ 5^0 sin 27rxi, 


^ = 2 1 ^ / 2 ^ (djW 

dt 3pdxi \ dxi J 3 p \dxi J ’ 


where a suitable nondimensionalization has been made, and go is a constant. 
In addition, p{6) and k{6) are, respectively, dimensionless forms of the vis¬ 
cosity and thermal conductivity and are functions of 9. Corresponding to 
the initial-boundary-value problem of the BGK model solved numerically, 
the above equations should be considered in the interval —1/4 < xi < 1/4 
with the Neumann conditions 


du2 

dxi 


= 0 , 


do 

dxi 


= 0 , 


at xi 



We should keep in mind that the compressible Navier-Stokes equations 
listed above hold only approximately because ui is not exactly zero in the 
numerical solution based on the BGK model. Now, we assume that p ~ 
const and du 2 /dt is negligibly small in the third equation. Integrating this 
equation with respect to xi and taking into account the boundary condition, 
we have du 2 ldx\ k, {g^pj^i^p) cos27rxi. We insert the expression of du 2 ldx\ 
in the forth equation and integrate it with respect to xi from xi = —1/4 to 
1/4 assuming that p ss const and p{6) ^ p{0av)- Then, we obtain 


(ft 


alp 

3(27r)2' 


Suppose that p{6av) = Cfj_9^y. Then, for the initial condition 6* = 1 at t = 0, 
we obtain 0^+'^ ps Cit -|- 1 with Ci = (1 -|- 5)gQp/12Tr‘^Cn, or for large t, 


Since d9/dt = 0{d9av/dt) = and p{du 2 /dxi)‘^ = 0(1/p{9av)) = 

(9(^-<5 /(i+'5 )), they are small in the energy equation. If we neglect these terms 
in the energy equation, we obtain dO/dxi ^ 0, i.e., 6 ss 9^^. Then, from the 
expression du 2 ldxi ss (gQp/2iTp) cos27rxi and the boundary condition, U 2 is 
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obtained as U 2 ~ {gopjAir'^/j,) sin27rxi. To summarize, we obtain 

9 ^ Cgt^, U 2 ~ sin 27rxi, (for t » 1), 

with new constants Cg and Cu- Because 5 = 1 for the BGK model, we 
have 6 ~ and U 2 ~ sin27rxi for large t. On the other hand, 

since 6 = 0.5 for hard-sphere molecules, 9 and U 2 behave as 0 ~ and 

U 2 ~ sin27rxi. Although this conclusion for hard-sphere molecules 

is consistent with the numerical result based on the Boltzmann equation, 
it does not coincide precisely with the numerical result based on the BGK 
model. It is natural because the argument is too sketchy. However, it pro¬ 
vides qualitative explanation for the slow increase of the temperature and 
its uniformity and for the slow decrease of the flow speed and the sinusoidal 
flow-velocity profile. More specifically, as the result of the temperature rise 
caused by the viscous heating, the viscosity and the thermal conductivity 
increase. The high thermal conductivity leads to a uniform 9. On the other 
hand, the high viscosity tends to prevent the external force from causing 
the gas flow, so that the flow speed decreases as the temperature increases. 
However, the decrease of the flow speed results in the decrease of the vis¬ 
cous heating. The rough estimate based on the compressible Navier-Stokes 
equations show that, as time proceeds, the amount of viscous heating de¬ 
creases, but the total amount of heat produced from the initial time increases 
indefinitely. Thus, the temperature continues to increase indefinitely. 


Conclusion 

We have discussed the numerous similarities between the steady prob¬ 
lem for the Boltzmann equation and for the Navier-Stokes-Fourier system. 
In particular, the presence of the viscous heating term in the tempera¬ 
ture equation depends on the scaling of the external divergence-free force 
field. However, we have observed a significant difference between the steady 
Navier-Stokes equation and the steady Boltzmann equation: while the steady 
Navier-Stokes equation with prescribed divergence-free external force always 
have a solution, either in the periodic setting or in the case of a bounded spa¬ 
tial domain with Dirichlet boundary condition, the steady Boltzmann equa¬ 
tion with nonzero divergence-free external force field cannot have a nonzero 
solution. We have proposed a physical explanation for this difference, based 
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on the long time behavior of the evolution Boltzmann equation with nonzero 
divergence-free external force field in the periodic setting. Our numerical 
computations based on the BGK model suggest that the temperature field 
in the gas increases indefinitely, so that the gas cannot reach a steady state. 

While we have presented our results on the hydrodynamic limit of the 
steady Boltzmann equation in the case of a hard sphere gas, the same results 
should remain true in the case of hard cutoff potentials (in the sense of Grad). 

Finally, our numerical simulations suggest the following problem, which 
we believe is open at the time of this writing. 

Problem: Gonsider the evolution Boltzmann equation set in the peri¬ 
odic box, with a prescribed nonzero divergence-free external force field (for 
instance f{xi,X 2 ,X 3 ) = (0,/o sin(27rxi), 0)). What is the asymptotic behav¬ 
ior of the temperature field averaged in the space variable in the long time 
limit? In particular, does there exist a > 0 such that the average temper¬ 
ature field is asymptotically equivalent to Ct" for some C > 0 as the time 
variable t tends to infinity? 

Note added after publication: After the present article was pub¬ 
lished by the Bulletin of the Institute of Mathematics of Academia Sinica, 
we became aware of the reference m This reference provides a complete 
proof of the asymptotic limit described in Theorem 1 at the formal level. 
The result established in |18) assumes that the source term in the Navier- 
Stokes-Fourier limiting system is small. 
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at Ecole polytechnique, and completed while the second author was visiting 
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Appendix: Properties of A and B 
Lemma 3. The tensor field A and the vector field B defined by the formulas 

A{v) :=u(8)u— B{v) := — 5)v 


satisfy the following properties. 
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(1) The orthogonality relations 

A _L Ker £ , i? _L Ker £ , A T B 
hold componentwise in L‘^ , Mdv). 

(2) There exists a unique tensor field A and a unique vector field B in 
L2(R3;(1 

+ \v\‘^)Mdv) such that 

CA = A and A _L Ker £ , 

CB = B and B _L Ker£ . 

(3) The tensor field A and the vector field B are of the form 

A{v) = a{\v\)A{v), B{v) = fi{\v\)B{v) 


for a.e. v £ R^. 

Statement (1) is Lemma 5.3 in [22], while statement (3) is Lemma 5.4 
in |22|. Statement (3) has been used systematically in the literature on the 
Boltzmann equation, referring to §7.31 in |12j . However, the discussion in 
[12| is incomplete. The key argument leading to the structure of A and 13 in 
statement (3) is the invariance of the linearized collision operator under the 
orthogonal group O3 (R). It seems that the first complete proof of statement 
(3) following this idea is in [T5|. Statement (2) follows from Hilbert’s lemma 
(the fact that the linearized collision operator £ satisfies the Fredholm al¬ 
ternative) . 

Lemma 4. The tensor fields A and A, and the vector fields B and B satisfy 
the following properties. 

(1) For each i,j = 1, 2,3, one has 

and 

(BiBj) = Sijl{\{\v\‘^ - 5)^|up/3(|u|)) = Kdij , 
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where 

= \{B-B) = \{B-CB)>^. 

(2) For eaeh i, j, k,l = 1, 2,3, one has 

{AijAki) = {SikSji + 6ii6jk - lhijSki)^{\v\*) 


and 


1/ :=^(|?;|^a(|u|)) 


where 

(3) For each i, j, k,l = 1, 2,3, one has 

{BiVjAj^i) ^ki{5i]^5ji F dii6jj^ 


and 


{BiUjAki) — c(^6ijf;6ji + diidjk , 


with 

c := ^{{B^v) : A) . 


Statement (2) is Lemma 3.4 in |22] (proved in Appendix 2 of [22]). See 
formulas (4.10) and (4.13b) in |5] for statement (1), which is proved by 
a similar (though simpler) argument as statement (2). (Notice the slight 
difference in normalization in the definitions of k in |5|and here). 


Proof of statement (3). By statement (3) in the previous lemma 

{BiVjAki) =(/l(|u|)|(|up - 5)viVjVkVi) 

- ^ki{P{\v\)l{\v\‘^ - 5 )\v\‘^ViVj) 
=(^(kl)5(I^P - 5)viVjVkVi) 

- hi{f3{\v\)l{\v\'^ - hfviVj) , 
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where the second equality follows form the fact that B _L Ker£. Then, 
according to Lemma 4.3 in [22], one has 

(/ 3 (| u |) 5 (| u |^ - 5)viVjVkVi) = X{6ij6ki + 6ik6ji + 6ii5jk) 

for some A G R, while, by statement (1) of the previous lemma, 

(/3(|u|)g(|up - 5)'^ViVj) = liBiBj) = Indij . 


Hence 


{BiVj T Sijf;6ji + ^u^jk) j 


and since 


(5A - 2K)dij = (BiVjAkk) = {BiVj trace(^)) = 0 , 

we conclude that A = |k, which immediatly implies the first equality in 
statement (3). The second equality is obtained in exactly the same manner. 
Contracting i, k and j, I, one finds the announced formula for c. □ 

Lemma 5. Let C be the tensor field defined by the formula 

C{v) := ® V ® V — 2>v ® I). 

Then, for each i, j, k,l = 1, 2,3, one has 

{Bi ^ Cjkl) — dikdjl ~\~ dilSjk^ . 


Proof. First 

{Bi ® Cjki) = \{BiVjVkVi) , 

since B T KerT. Thus 

\{BiVjVkVi) = - 5)viVjVkVi) 

hidijSki T dikdji T diiSjk) 

for some /r G R, by Lemma 4.3 in |22] . By contraction on the indices k, I 
bfidij = liBiUjlvl"^) = {BiBj) = ^ij 
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where the first equality follows from the orthogonality relation B _L Ker£ 
by statement (1) in LemmaHence fj, = ^, and this immediatly implies 
the desired identity. □ 
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